Glassy dynamics, metastability limit and crystal growth in a lattice spin model 
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We introduce a lattice spin model where frustration is due to multibody interactions rather than 
quenched disorder in the Hamiltonian. The system has a crystalline ground state and below the 
melting temperature displays a dynamic behaviour typical of fragile glasses. However, the super- 
cooled phase loses stability at an effective spinodal temperature, and thanks to this the Kauzmann 
paradox is resolved. Below the spinodal the system enters an off-equilibrium regime corresponding to 
fast crystal nucleation followed by slow activated crystal growth. In this phase and in a time region 
which is longer the lower the temperature we observe a violation of the fluctuation-dissipation theo- 
rem analogous to structural glasses. Moreover, we show that in this system there is no qualitative 
difference between a locally stable glassy configuration and a highly disordered polycrystal. 



When cooled fast enough, most liquids can be brought 
below their melting temperature T rn while avoiding crys- 
tallization. In this supercooled phase the viscosity in- 
creases with decreasing temperature, until at the glass 
transition T g relaxation time becomes too long and equi- 
libration cannot be achieved within experimentally ac- 
cessible times. In 1948, W. Kauzmann [[j] noted that if 
the entropy of a supercooled liquid is extrapolated below 
T g , it becomes equal to the crystal entropy at a temper- 
ature T s > 0, and even smaller than zero if extrapolated 
further. This entropy crisis is never actually observed, 
because the glass transition intervenes before. However, 
Kauzmann found it paradoxical that it was just a kinetic 
phenomenon (the glass transition) that saved the liquid 
from a thermodynamic nonsense. 

In the context of the Adams-Gibbs theory j| , the en- 
tropy crisis has however an interpretation: the entropy 
difference between crystal and liquid is related to the con- 
figurational entropy S, that is the entropic contribution 
due to the presence of an exponentially high number of 
different glassy minima. The vanishing of £ at T s sig- 
nals a thermodynamic transition to a new phase, char- 
acterized by a sub-exponential number of glassy states, 
separated by infinite free-energy barriers. This picture is 
exact for some mean- field spin-glass systems ||, and it 
may be the correct resolution of the Kauzmann paradox 
even for real structural glasses. 

Despite analytic and numerical work supporting the 
entropy crisis scenario Q, there is another way to avoid 
the Kauzmann paradox, which, interestingly enough, was 
proposed by Kauzmann himself jjj. He rejected the idea 
of a thermodynamic glassy phase, and of a transition at 
T s . What Kauzmann hypothesized is the existence of a 
metastability limit of the supercooled liquid phase, be- 
low which crystal nucleation becomes faster than liquid 
equilibration. More precisely, he defined an effective spi- 
nodal temperature T sp > T s below which "the free energy 
barrier to crystal nucleation becomes reduced to the same 
height as the barrier to simpler motions". Below T sp the 
supercooled liquid is operationally meaningless and thus 



the paradox is avoided. However, the metastability limit 
may prove impossible to observe experimentally if the 
equilibration time at T sp is much larger than the experi- 
mental time, that is if T sp < T g . 

In the Kauzmann scenario, the off-equilibrium phase 
below T sp basically consists in a very slow crystal domain 
growth. This suggests two criteria to detect whether a 
system has such a metastability limit or not, even when 
T sp is experimentally inaccessible. First, we may think 
that there is a qualitative difference between a disordered 
glassy configuration, obtained by quenching a liquid, and 
a polycrystalline configuration, however rich in defects 
this is, and however slow crystal growth may be. Sec- 
ondly, we know that the off-equilibrium dynamics of sim- 
ple domain growth (as in the Ising model) can be dis- 
tinguished from glassy dynamics (as observed in struc- 
tural glasses) by a different violation of the fluctuation- 
dissipation theorem (FDT) j| ||. 

In this Letter we present a model where the Kauzmann 
paradox is avoided by a metastability limit. The model is 
thus a good test for the two criteria discussed above. We 
shall find that neither criterion is sharp enough to dis- 
criminate such a system from a typical structural glass. 
More precisely, if our experimental time were not long 
enough to explicitly observe the loss of stability of the 
liquid at T sp , it would be impossible to distinguish the 
present system from an ordinary fragile glass. The rea- 
son is that below T sp crystal nucleation is fast, but crys- 
tal growth becomes very slow, with many crystal droplets 
trying to expand in a liquid background |7j. In such a 
situation distinguishing between a truly disordered glass 
and a mixture of tiny mismatched crystallites becomes 
very hard, and FDT violation is nontrivial. 

The aim of our study is twofold. First, we want to show 
that Kauzmann's resolution of the Kauzmann paradox is 
valid at least in one simple system, and that it is not nec- 
essarily in conflict with glassy phenomenology. Second, 
we hope this example will help to develop some more 
strict criteria to distinguish systems with a thermody- 
namic transition at T s , from those where the Kauzmann 
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scenario holds. The Hamiltonian of our model is, 

N 

H = J2(l + m)fi, fi = sT r sfs?s», (1) 

i=l 

where W is for west, S for south, etc. The spins Sj = ±1 
belong to a two-dimensional square lattice of linear size 
L. We perform single-spin- flip Monte Carlo simulations 
in square lattices with L = 100 and L = 500 ||. 

The disordered version of this model was first intro- 
duced in H to describe an ensemble of coupled two level 
systems, whence its name, CTLS. There has been re- 
cently much interest in lattice models of this sort, where 
multibody interactions in the Hamiltonian ensure frus- 
tration even without quenched disorder [jl(| O, |l2], |l|. 

The CTLS has a crystalline ground state obtained by 
covering the lattice with the following non-overlapping 5- 
spins elements: Sj = —1, — sf = sf = sf = +1. The 
ground state energy density is ecs = —1-6. The unit 
cell size is 5 x 5, and this together with the symmetry 
x — > T x ,y — ► ±2/, gives a ground state degeneracy of 
50. The finite degeneracy of the crystal is a key feature 
of the CTLS compared to the plaquette model of pi] , 
since it will allow us to directly measure the amount of 
crystalline order in the system. 

Model (||) has a first order (melting) transition. To 
locate the melting temperature T m we compute the free 
energy / for the crystal (CR) and liquid (LQ) and ask 
that /cR(Tm) = fhQ(T m ). The free energy is obtained 
integrating the energy (/3 = 1/T): pf(fi) = #,/(#,) + 
fpd/3' e(/3'), taking /3 = for the liquid, and /3 = oo for 
the crystal. The equilibrium energies are well fitted by, 

e LQ (T) = -1.8 tanh(l.l/T) , (2) 
ecR(T) = e GS + 5.5 x lO" 5 T 13 . 

Using these relations we find T m = 1.30. Unless cool- 
ing is exceedingly slow, crystallization is not attained at 
T m , and the liquid can be kept at equilibrium in its su- 
percooled phase T < T m . Extrapolation of clq{T) and 
/lq (T) gives the temperature T s where the entropy of the 
supercooled liquid equalizes that of the crystal, namely 
the Kauzmann paradox temperature. We find T s = 0.91. 

The equilibrium dynamics of the CTLS in the super- 
cooled phase can be studied by measuring the normalized 
spin-spin correlation function, C n (t,t w ) = [(s(t w )s(t)) — 
(s) 2 ]/[l - (s) 2 ], with t > t w . In equilibrium C n (t,t w ) — 
C n (t — tw), and we find that the correlation can be fit- 
ted to a stretched exponential, C n = exp [— (i — itu/r)' 3 ] . 
In Fig. [l] we plot the relaxation time r as a function of 
T, together with a power law fit r = A/(T — T^) 7 , with 
T c = 1.06, 7 = 2.29. The accuracy of the fit suggests 
that the CTLS is a fragile system, a fact supported by 
the Angell plot jl4| in the inset of Fig. |l|, which com- 
pares the relaxation times of the CTLS and a strong 
system, the two-dimensional plaquette model (2d-PQ) 




FIG. 1: Relaxation time as a function of the temperature. Full 
line: power law fit. Inset: fragility plot for a comparison of 
the CTLS with the model studied in [[L2| . T g is an operational 
'glass transition', defined by r(T 9 ) = 1000. L = 500. 

studied in p2| . Data are fitted to a Vogel-Fulcher form, 
t = t exp[A/(T - T )}. In the CTLS we find T = 0.76 
for T G [1.2 : 1.7], and T = 0.90 for T G [1.2 : 1.4]. 

Below about T w 1.2 the relaxation time cannot be 
measured, because crystal nucleation starts before the li- 
quid equilibrates. This can be seen in Fig. where we 
plot the energy density vs. time at four different tem- 
peratures below T m , with random initial condition. At 
T = 1.26, the system relaxes in the supercooled liquid 
and remains in this phase up to our experimental time, 
2 x 10 6 MCS. On the other hand, for T = 1.23, after ini- 
tial equilibration in the liquid phase, the system makes 
a sharp transition to the crystal: on average crystal nu- 
cleation starts at about 10 5 MCS, while crystal growth 
is completed in 10 6 MCS. At T = 1.18, however, nucle- 
ation starts roughly after 10 4 MCS, while complete crys- 
tallization is achieved in more than 10 6 MCS. At T = 1, 
the liquid plateau disappears and crystallization is still 
incomplete after the largest time. Thus at lower temper- 
atures nucleation is faster, but crystal growth is slower. 

To pinpoint the liquid metastability limit we use stan- 
dard nucleation theory |1^| . A necessary condition for 
the existence of the supercooled liquid is that the crys- 
tal nucleation time r nuc is much longer than the liquid 
equilibration time r eq . We can estimate t„„ c (T) close to 
a reference temperature T*, assuming that the surface 
tension is constant. We have, 

r nuc (T) = exp i ^|y log[r„ uc (T*)]} , (3) 

where 6f(T) = /lq (T) - f CR (T) is the bulk free energy 
difference between supercooled liquid and crystal. Wc 
have chosen T* = 1.23. To estimate the liquid equili- 
bration time we note that the correlation C{t) drops to 
zero in about 20 relaxation times, and thus r eq = 20 r. 
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FIG. 2: Left: e vs. t at different temperatures, L — 100. 
Error bars are only showed when larger than symbols size. 
The values of m and £ in the legend are measured at the 
longest time for each temperature (see text). Right: crystal 
nucleation time and liquid equilibration time vs. temperature. 



By imposing T nuc (T sp ) = T eq (T sp ), we obtain the effec- 
tive spinodal temperature T sp marking the metastability 
limit. Fig. | (right) shows that T sp = 1.22 jjTj. We have 
T sp > T s , and the Kauzmann paradox is thus avoided. 

Below the metastability limit T sp the only equilibrium 
phase is the crystal. However, a long-lived off-equilibrium 
glassy phase can still be formed: if we cool the system, 
it eventually remains stuck in an off-equilibrium state 
whose asymptotic energy is lower the slower the cooling 
rate r (Fig. ||). All the configurations reached at T = 
are stable minima, indicating that activation is needed to 
grow the crystal. For very fast coolings these configura- 
tions are completely disordered, while at the slowest cool- 
ing rates they correspond to highly ordered polycrystals. 
Is there a qualitative difference among these asymptotic 
states ? In other words, is it possible to sharply separate 
a bona fide glass from the polycrystal ? 

To answer this question we introduce two direct mea- 
sures of crystalline order. The first one is the normalized 
domain size £, obtained from the Fourier transform G(k) 
of G(r) — (siS i+r ). G{k) has a peak at k = 2n/5, and 
we define £ as the inverse of the peak width, normal- 
ized by L. The second is the normalized crystal mass 
to, i.e. the total fraction of crystallized spins (defined 
as the number of spins down surrounded by 8 spins up, 
multiplied by 5) divided by L 2 [jl6|]. Using £ and m we 
see that crystallization has been achieved at the lowest 
three temperatures of Fig. ||. If we now plot the energy of 
the T — asymptotic configurations in a cooling exper- 
iment vs their crystalline mass and domain size, (Fig. 
inset), we find a continuous spectrum of states. Thus, in 
the CTLS there is structural continuity between highly 
disordered glassy minima and strongly ordered polycrys- 
talline minima. The answer to the question above is 
therefore negative. Yet, had our slowest cooling rate been 




T 



FIG. 3: e as a function of T, at various cooling rates r. Full 
line: equilibrium liquid energy, eq. (0). Inset: e vs m and £ 
in the T = asymptotic states (r £ [10~ 6 , 10" 3 ]). L = 100. 

r = 2 x 10~ 4 we would only observe a disordered glass. 

If not by structural difference, we might expect the 
dynamic behaviour of our "glass" to betray its nature 
of crystal-growth phase, by displaying a characteristic 
coarsening dynamics. In coarsening, excess energy over 
the ground state is concentrated in the interfaces among 
domains. This gives e(t) - e G s oc £(*) d ~7£(*) d = 1/f (t), 
for L£ 3> a, where a is the interfacial width. We do in 
fact find such a regime (Fig. ||) for £ > 0.2. However, due 
to slow activated dynamics, the early regime, L£ ~ a, can 
last for long, and in fact it is the only one that can be 
observed at the lowest temperatures. This regime, which 
we call bubbling, is characterized by a rapid increase of 
to at roughly zero £ (Fig- inset), meaning that fast 
nucleation leads to rapid formation of many tiny crystal 
droplets. Only at longer times, when most of the system 
has crystallized, the domains grow at expense of each 
other and proper coarsening starts. For large systems we 
expect m ~ 1 — a/L£, such that bubbling and coarsen- 
ing regimes become well separated. Summarizing, at low 
temperatures, or for short experimental times, the coars- 
ening regime is inaccessible, in the same way as, for fast 
cooling, polycrystals are not observed. 

This behaviour has important consequences on the 
FDT violation pattern. The integrated response is 
given by x{t,tw) = J^dt' R(t,t'), where R(t,t') = 
6(s(t))/6h(t') and h is the conjugate field to the spin. We 
can make a parametric plot of x(t, t w ) vs the correlation 
C{t,t w ) = (s(t)s(t w )). At equilibrium FDT holds and 
X = — C), while a departure from FDT is normally 
observed for late times if the system is out of equilibrium. 
This FDT violation is the standard tool to discriminate 
genuine glassy behaviour from simple domain growth. 
The relation x — 0X(l — C) can be used to express the 
late time FDT violation in off-equilibrium systems. For 
simple domain growth X = 0, while in structural glasses 
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FIG. 4: Excess energy e — eGS as a function of £. Full 
line: l/£ fit of the data. Inset: m(t) as a function of {(t), 
parametrically in t, after a quench at various temperatures. 
All runs are 2 • 10 6 MCS long. L = 100. 




FIG. 5: Parametric plot of x(t>t-ui) vs C(t,t w ), at various val- 
ues of the waiting time t w . Full lines represent the equilibrium 
relation \ = — C). For T = 0.66 the broken line is a linear 
fit for C < 0.8. L = 500. 



it is observed < X < 1 @, f| 18 . We do expect a 
trivial FDT violation (X = 0) in the coarsening regime, 
but the situation may be more complicated in the earlier 
bubbling regime. From Fig. |^ we see that at T = 0.66 the 
system remains in the bubbling regime up to the largest 
time. At the same temperature we find an FDT violation 
very similar to structural glasses (Fig. [|): for C smaller 
than a given breaking point the slope changes, giving rise 
to a constant nonzero value of X. Moreover, the pattern 
basically does not change within two orders of magnitude 
in the waiting time t w . At T = 0.66 we can thus define an 
effective temperature T e ft = T/X = 1.32. On the other 
hand, for T = 1 as t w increases we see a crossover from a 
nontrivial FDT violation (X ^ 0), to a seemingly trivial 
pattern (X ~ 0). This is consistent with the fact that at 
T = 1 the coarsening regime is reached within our exper- 



imental time (Fig. ||). Yet, once again, had our longest 
experimental time been too short to enter the coarsening 
regime, we would only observe a glassy FDT violation. 

In this Letter we presented a model where the metasta- 
bility limit can be either observed directly as a loss of sta- 
bility of the supercooled liquid, or indirectly, thanks to 
the formation of polycrystalline asymptotic states in slow 
cooling experiments, and thanks to trivial FDT violation 
in the late coarsening regime. However, in order to be 
effective the indirect methods need as long an experimen- 
tal time as the direct method. For shorter time scales the 
present model is compatible with ordinary fragile glasses. 
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